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Background 


• Boundary Layer Ingestion (BLI) Propulsion has the potential 
for significant reduction in Aircraft Fuel Burn (5-10%) 

• Previous studies referenced in 2011 FAP presentation by 
Razvan Florea: 

Bangert, etal., NASA-CR-3743 (1983) 

Daggett, etal., NASA-CR-2003-21 2670 
Berrier, NASA-TP-2005-2 13766 
Campbell, AIAA 2005-0459 
Kawai, etal., NASA-CR-2006-2 14534 16% 

Carter, AIAA JOA 2006, Vol 43, No. 5 
Plas, MIT PhD Thesis 2006 
Plas, et at., AIAA 2007-450 10% 

Kawai, NASA-CR-2008-215141 

Burn Reduction 

Nikol, NASA-TM-2008-215112 

5% 

Drela, AIAA 2009-3762 
Nikol, McCuller, AIAA 2009-931 
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Boeing Boeing MIT NASA 

(Daggett, et al. f 2003) (Kawai, et al., 2006) (Plas, et al., 2007) (Nickol, et al., 2009) 


Technical Challenges 

• The potential benefits of Boundary Layer Ingestion (BLI) 
Propulsion can be diminished by considerations of 

- Inlet total pressure loss 

- Fan efficiency reduction 

- Fan stall margin reduction 

- Fan aeromechanics 

(dynamic stresses and flutter stability) 

V / 
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Optimization-Based Parametric Inlet Design 


• Inlet excess pressure loss reduced ~4-5x relative to original Inlet A starting point 

• Dominant distortion harmonic amplitudes reduced ~ 30-50 % relative to original 
^ Inlet A starting point 

5 
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Fan Efficiency with Distortion-Optimized Inlet 


Inlet significantly 
improves fan interaction 
with incoming distortion 


Fan Efficiency Reduction (%) 



Distortion Index (ma< AP C / P, %) 

Fundamental Aeronautics Program 
Subsonic Fixed Wing Project 


Excursions in Fan Blade Leading Edge 
Relative Incidence from Clean Inflow 


Inlet A P3 Inlet 
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Inlet enables fan to 
meet performance 
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Fan CFD Analysis - TURBO Code 


Implicit, finite-volume solver 

Reynolds-Averaged Navier Stokes equations 

Structured multi-block code 

Multi blade-row code 

k-epsilon turbulence model 

Inlet distortion boundary condition 

Throttle exit boundary condition 



Dynamic grid deformation for blade vibration 

Prescribed harmonic blade vibrations with energy method 
to evaluate flutter stability 
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Fan Computational Domain 


Analysis of an Aero Design Iteration (not the Final Design) 



Fan rotor 

case 

blade 



Exit 

boundary 


Core-bypass 

splitter 


Exit 

boundary 


H, O, and C blocks 
of mesh generated by UTRC 
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Total Pressure Ratio 


Fan Performance - Clean Inflow 



• TURBO code (RANS solver) used with radial inlet profile 
of total pressure, total temperature, and flow angles 

• Speedline traversed by setting exit throttle condition and 
converging flow solutions 
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Inlet Flowfield Provides Distortion Pattern 


• Inlet flow computations were performed at UTRC for an 
inlet design iteration (not final design) and the flowfield 
results were provided to NASA 
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Fan Computation with Inlet Distortion 




Inlet distortion is prescribed as boundary condition at inlet 
boundary of the fan computational domain ( 18 -blade fan 
rotor and splitter) 
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Total Pressure Ratio 


Periodicity of Flowfield Around the Rotor 


• Total pressure ratio for various blade passages 



Fundamental Aeronautics Program 
Subsonic Fixed Wing Project 


Time Step Counter 


Total Pressure Ratio 


Periodicity of Flowfield Around the Rotor 


• Total pressure ratio for various blade passages 

1.45 
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1.25 
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passage 10 



Inlet Distortion causes variations in mass flow 
rate and pressure ratio around the fan rotor 
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Aeroelastic Formulation 



• Blade structural dynamics modal equations with 
aerodynamic load 

+[*]{<?} = { ad} 

aerodynamic load vector - 

Modal Force 



\ 

Modal Force computation 
requires unsteady pressure 
and modal displacements 





Forced Response 
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Structural Dynamics Model & Results 


Blade structural model created based on aero 
design iteration (structural design is in progress) 

• 8-node brick elements 

• 9,782 elements, 15,096 nodes 

• 222 nodes at the root constrained 



mode 1 mode 2 mode 3 mode 4 

63.5 Hz 156.6 Hz 224.8 Hz 346.6 Hz 
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Modal Force 


Time history over one rotor revolution 
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Time Step Counter 


Magnitude, lb 


Modal Force 


Fourier components 


700 


600 


500 


400 


300 


200 


100 


mode 1 
mode 2 
mode 3 
mode 4 



1 23456789 10 


Fundamental Aeronautics Program 
Subsonic Fixed Wing Project 


Harmonic Number 


Frequency (Hz) 


Campbell Diagram 
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Magnitude, lb 


Modal Force 


Fourier components 
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Harmonic Number 


generalized displacement 


Forced Response - 

Vibrat i on A mplitude and Dynamic Stresses 



• Dynamic stresses are required to determine fatigue 
characteristics (Goodman diagram) 


{^nr} [[^n] 

dynamic stress 



l 


{ AD nr | for n th mode, r th harmonic 


where s n is the modal stress 



harmonic or 
engine order 

vibration amplitude 
(inch) at tip t.e. 

dynamic stress 
amplitude (psi) 

1 

5.5x10 2 

273 

2 

3.0x10- 2 

290 

3 

1.9x10-2 

666 

4 

3.1 x 10- 3 

308 

5 

2.6 xIO- 3 

169 

6 

2.7 xIO- 4 

33 

7 

7.0x1 O' 4 

427 

8 

6.0x10- 5 
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Flow Chart for Flutter Stability Computation 



Aerodynamic damping computation using TURBO-AE 



Configuration 
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Flutter Stability with Clean Inflow 



Design operating speed, mode 1, 0 nodal diameter pattern 
(all blades in-phase), 18 blade passages (full rotor) 
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Flutter Stability with Clean Inflow 


• Design operating speed, 18 blade passages (full rotor) 

• Phase angle of vibration = 360 * Nodal Diameter / 1 8 
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Flutter Stability with Clean Inflow 



• Design operating speed, 18 blade passages (full rotor) 

• Phase angle of vibration = 360 * Nodal Diameter / 1 8 
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Aero Damping 


Flutter Stability with Clean Inflow 


• Design operating speed, 18 blade passages (full rotor) 

• Phase angle of vibration = 360 * Nodal Diameter / 1 8 
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Fast-Running Aeroelastic Analysis 



• Harmonic Balance CFD Method 

• Fourier series expansion 

N 

U(x,t)~ '2 / U n (x)e i " m 

n=-N 

substituted into governing equations and solved for each 
harmonic component [Hall, 2000] 

• Lax-Wendroff method 

• 2 nd and 4 th order smoothing for stability 

• Non-reflecting boundary conditions 

• Spalart-Allmaras turbulence model 

• Eigenvalue analysis to calculate aerodynamic damping 
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Aerodynamic Damping 


Harmonic Balance Results - Clean Inflow 
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Flutter Stability with Distorted Inflow 


Various Approaches 

• Circumferentially average the distorted inflow to obtain an 
equivalent radial profile; use work-per-cycle analysis 

• Select a portion of the inlet distortion to represent a “worst- 
case” inflow condition that is used at all circumferential 
locations; use work-per-cycle analysis 

• Prescribe blade vibrations and distorted inflow; use work- 
per-cycle analysis; average the results over all blades, and 
over multiple blade vibration cycles 

• Use tightly-coupled aeroelastic analysis with distorted 
inflow; blade vibrations are determined as part of the 
computations; post-process time history to estimate 
average damping over all blades and multiple vibration 
cycles 
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Flutter Stability with Distorted Inflow 



Current Preferred Approach 

• Prescribe blade vibrations and distorted inflow 

• Use work-per-cycle analysis 

• Average the results over all blades, and over multiple blade 
vibration cycles 



cycle J surface 


-p.d A • ( )dt 



dt 




Unsteady pressure includes effect of 



isolate this component to 
assess flutter stability 


2) blade vibration i=> 


v. 
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Flutter Stability with Distorted Inflow 



Design operating speed, mode 1, 0 nodal diameter pattern 
(all blades in-phase), 18 blade passages (full rotor) 
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Vibration Cycle Number 
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Summary 


• Created structural model based on aero design iteration 
and computed structural dynamics characteristics 

• Performed aeromechanical analysis of design iteration 

• Performed fan flutter analysis with clean inflow at design 
speed - no flutter encountered at conditions analyzed; 
additional work needed at part-speed conditions 

• Performed distorted inflow analysis for forced response 
vibrations to determine dynamic stress at design speed - 
additional work needed at on-resonance conditions near 
design speed 

• Performed initial analysis with blade vibrations and 
distorted inflow to estimate flutter stability - additional 
flutter analyses needed for other vibration modes and 
operating conditions 
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Future Work 


• Extend computational domain to include fan exit guide 
vanes in the unsteady aerodynamics analysis 

• Perform aeromechanical analysis on updated fan stage 
design with non-axi-symmetric exit guide vanes 

• Perform aeromechanical analysis on final inlet-fan design 
to ensure safe wind-tunnel test 

• Develop tightly-coupled aeroelastic analysis capability in 
TURBO for more detailed analysis of blade vibrations with 
distorted inflow 

• Develop inlet-fan coupled aeroelastic analysis capability 
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